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We analyze numerically the performance of the near-optimal quadratic dynamical decoupling (QDD) single- 
qubit decoherence errors suppression method [J. West et al, Phys. Rev. Lett. 104, 130501 (2010)]. The QDD 
sequence is formed by nesting two optimal Uhrig dynamical decoupling sequences for two orthogonal axes, 
comprising Ni and N2 pulses, respectively. Varying these numbers, we study the decoherence suppression 
properties of QDD directly by isolating the errors associated with each system basis operator present in the 
system-bath interaction Hamiltonian. Each individual error scales with the lowest order of the Dyson series, 
therefore immediately yielding the order of decoherence suppression. We show that the error suppression prop- 
erties of QDD are dependent upon the parities of Ni and iV2, and near-optimal performance is achieved for 
general single-qubit interactions when iVi = N2. 
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I. INTRODUCTION 

In recent years, there have been promising advances to- 
wards the usage of quantum systems to perform quantum in- 
formation processing (QIP) 0"). However, for these systems 
to be utilized efficiently, it is necessary to preserve and store 
information for a sufficient amount of time so that computa- 
tions can be implemented. Unfortunately, quantum systems 
are generally hindered in their ability to perform such tasks 
effectively due to unwanted interactions between the system 
and its enviroment, which result in decoherence J5) . 

Dynamical Decoupling (DD) is a strategy that can be used 
to suppress decoherence and effectively remove unwanted 
system-environment interactions for a period of time such 
that quantum states can be preserved with a low probabil- 
ity of error Q. As originally conceived, DD schemes are 
characterized by the application of short control pulses to the 
system such that the overall time evolution provided by the 
pulses selectively average out system-environment interac- 
tions, thereby suppressing decoherence EHT]. One advantage 
of DD is that it is an open-loop control method, which does 
not require any measurements or feedback, unlike quantum 
error correction |8 1. Nor does DD require any specific knowl- 
edge of the environment other than it being non-Markovian, 
unlike optimal control methods designed to suppress decoher- 
ence ll9l[T0l. 

Early DD schemes were designed to remove unwanted 
system-bath interactions to a given, low order in time- 
dependent perturbation theory [3 |. Concatenated DD (CDD) 
was the first explicit scheme capable of removing such inter- 
actions to an arbitrary order [11]. CDD accomplishes this 
via a recursive construction in which each successive level 
removes another order in time-dependent perturbation the- 
ory. The advantages of CDD over standard periodic pulse 
sequences have been extensively studied analytically lfT2l [T3l 
and numerically lfl2l [Hl - flTl . and confirmed in a number of 



recent experimental studies 1 18-22]. However, assuming that 
pulse intervals can be made arbitrarily short, the number of 
pulses required to achieve arbitrary order suppression grows 
exponentially with the order in CDD. When the finiteness of 
pulse intervals is accounted for there is an optimal level of 
concatenation and correspondingly a highest attainable per- 
turbation theory order for removal of unwanted interactions 

nuns ed. 

In contrast to CDD, Uhrig DD (UDD) is characterized by 
the use of unequal pulse intervals, or free evolution periods 
||23ll . By applying control pulses at 
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where j = 1,2,...,7Y+1, a UDD sequence of total dura- 
tion T yields Ath order decoupling for single-axis system- 
bath coupling l24l . The number of pulses comprising a UDD 
sequence is A if A is even, or A + 1 if A is odd. General- 
izations of UDD for generic system-environment interactions 
include Concatenated UDD (CUDD) (25], Quadratic Dynam- 
ical Decoupling (QDD) [26] and Nested Uhrig Dynamical 
Decoupling (NUDD) E3 (see also |28)). CUDD removes 
the restriction of single-axis decoupling and suppresses gen- 
eral (three-axis) decoherence errors on a qubit, but still suf- 
fers from the exponential cost of CDD. On the other hand, 
QDD also suppresses general decoherence errors on a qubit, 
but does so without the exponential cost of CDD by nesting 
two UDD sequences for two orthogonal axes. In fact, QDD 
is a near-optimal scheme for single-qubit decoupling from an 
arbitrary bath, requiring only A 2 pulses for Ath order decou- 
pling if A is even, or (A + 1) 2 pulses if A is odd. The NUDD 
sequence is built on the same nesting idea, but removes the 
QDD restriction of single-qubit decoupling. NUDD applies 
to arbitrary multi-level systems coupled to arbitrary baths, as 
recently proved in Ref. ll29ll . 

In this work we focus on the decoherence errors suppres- 
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sion capabilities of QDD. This protocol requires two nested 
sequences each containing Nj pulses if Nj is even or Nj + 1 
pulses if Nj is odd, j — 1,2. We call N% and N 2 the inner and 
outer sequence order, respectively. While the original QDD 
paper |26| noted that sequences with Ni ^ N 2 are possible 
and could be advantageous when a particular axis is dominant, 
only the case Ni — N 2 was analyzed. Here we numerically 
study QDD for N% ^ N 2 . We provide a complete numerical 
elucidation of the performance of QDD as a function of iVi 
and N 2 . 

The structure of this paper is as follows. In Section [TT] we 
provide a brief synopsis of our results. In Section III we sum- 
marize the QDD protocol. In Section IV we introduce the 



relevant error measures utilized in this paper and provide an 
analysis of the expected scaling of these measures with the in- 
ner and outer sequence orders. Section |V]is devoted to our 
numerical results. In it we discuss the scaling of the single- 
axis errors and their time-dependence, as well as the scaling of 



a distance measure for the entire QDD sequence. Section VI 
presents our conclusions. 



holds for the outer sequence as well (i.e., the a z interaction, 
except in the case of N 2 > 2(Ni + 1) with Ni odd, where 
the inner sequence hinders the ability of the outer sequence 
to suppress decoherence to the expected order. We find that 
in general, UDD efficiency for general single-qubit errors is 
achieved when N\ and N 2 are both even. 

We introduce another perspective on the QDD sequence, by 
studying the time-dependence of the single-axis errors. We 
find that as the sequence progresses, these errors oscillate be- 
tween values that are near to their final minimum, and much 
higher values. 



III. QDD PROTOCOL 

UDD suppresses single-qubit dephasing or longitudinal re- 
laxation errors separately, for a general environment 11301 [3D . 
The extension of UDD to QDD improves on this by handling 
general single-qubit decoherence, in particular both dephas- 
ing and relaxation simultaneously. Therefore, in our analysis 
of QDD the time-independent Hamiltonian 



n. SYNOPSIS OF RESULTS 

We characterize QDD performance with respect to the or- 
der of error suppression given by the overall fidelity loss of 
the qubit (system) state. We show that the order of overall er- 
ror suppression is dictated by the lowest of the inner or outer 
sequence orders, namely D ~ ( j T ) I * 1 ui{iVi,jV2}+i > w here £) j s 
a measure of the overall error, J is the strength of the system- 
bath coupling, and t is the smallest pulse interval. 

To arrive at these results and gain more insight we first char- 
acterize QDD performance with respect to the order of error 
suppression given by the order of error suppression along each 
axis of the qubit Bloch sphere, referred to as the single-axis 
error. We isolate the single-axis errors by projecting the to- 
tal evolution operator into the three directions defined by the 
Pauli basis. The error suppression properties are then distin- 
guished with respect to the scaling of the error as a function 
of the minimum pulse interval for various inner and outer se- 
quence orders. Since this scaling is dominated by the first non- 
zero term of the Dyson series, the order of error suppression 
for each single-axis error can be characterized with respect to 
inner and outer sequence order. Suppression of the x-axis er- 
ror (or "X-type error") is determined by Ni, the y-axis error 
by the parity of N% and N 2 , and the z-axis error by N 2 except 
in the case when Ni is odd [see Eqs. ([3~5]>-(|3"7]i]. Our results 
will show that if Ni is odd, there is a constraint on the sup- 
pression of the z-axis error depending on the value of Ni with 
respect to N 2 . 

Parity effects were anticipated in Ref. [27], where the ex- 
pected performance of QDD was proved for sequences with 
even Ni (the proof was recently completed in Ref. ||29l ). We 
show that parity effects are absent in QDD only for the inter- 
action that anti-commutes solely with the decoupling operator 
comprising the inner sequence. This interaction (a x ) is sup- 
pressed with UDD efficiency, i.e., A^tii order error suppres- 
sion for a pulse sequence comprising iVi pulses. This result 



H = H B + H SBl (2) 
H B = I®B T , (3) 
H S b = o x ® B x + a y ® B y + a z ®B Z (4) 

is employed to describe general system-environment interac- 
tions for the single-qubit system. The system operators are 
the standard Pauli matrices, while the bounded operators, B^, 
[i € {/, x, y, z}, characterize a generic environment. The op- 
erator Bi encompasses the pure bath dynamics, so that H B is 
the "pure-bath" Hamiltonian, while Hsb is the system-bath 
interaction Hamiltonian. 

In building a QDD sequence to compensate for the in- 
teractions present in the above Hamiltonian, it is useful to 
first choose a so-called Mutually Orthogonal Operator Set 
(MOOS) l27l . This set comprises mutually anti -commuting 
or commuting operators needed to perform DD in the gen- 
eral NUDD scheme. Each of these operators is required to 
anti-commute with some portion of the interaction Hamilto- 
nian. The anti-commutation condition between an element of 
the MOOS and the interaction Hamiltonian is essential for de- 
coherence suppression. One might expect a UDD sequence 
composed of an element of the MOOS to suppress the cor- 
responding anti-commuting interaction with UDD efficiency. 
As we shall show, this is not the case, except for the inner 
sequence. 

In the case of a single-qubit system, NUDD reduces to 
QDD and the MOOS requires only two operators. There is 
some flexibility in choosing the MOOS for Eq. but with- 
out loss of generality we pick fl = {X, Z}. The operators 
X = <j x ® I and Z = a z ® I (we drop global phase fac- 
tors everywhere) represent ideal zero-width it -rotations about 
their respective axes of the qubit subspace, and do not affect 
the bath. The QDD sequence is now readily constructed as 

ma 

U^' N2) =X N *zW{s Na+1 T)-XZW{s 1 T), (5) 
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where 

ZW( T ) = Z n 'U(s Ni+1 t) ■ ■ ■ ZU(s 2 t)ZU{ Si t). (6) 

The free evolution dynamics between successive pulses, 
U(t) = er %Ht , is governed by Eq. (|2j, and the free evolu- 
tion time durations are given in terms of the normalized UDD 
intervals, 



ti — to 



(7) 



with tj specified by Eq. (JTJ, and r = t\ — to- The high ef- 
ficiency of UDD in suppressing decoherence, and therefore 
QDD, arises from the choice of the relative free evolution time 
durations {sj}. 

The total normalized time of an TV-pulse UDD sequence is 
given by, 



N+l 

SfW EE £ Sj = 
3=1 



tN+1 2 
= CSC 



2N + 2 J ' 



(8) 



so that the total physical time is 



T (N) _ S (N), 



(9) 



Therefore the total normalized time of a QDD sequence 
with Ni inner and N 2 outer pulses is given by, 

N 2 + l 

3=1 



CSC 



esc 



(10) 



2Ni+2J \2N 2 + 2 
so that the total physical time is T^ Nl ' N ^ = S( Nl ' N *>T. 

IV. QDD PERFORMANCE MEASURES 

A. QDD analysis 

Our goal here is to understand the properties of QDD error 
suppression for general N% and N 2 and for a wide range of 
parameters. We do so by isolating the errors proportional to 
each system basis operator, a^, // = {x, y, z}, i.e., the single- 
axis errors. In this manner the order of error suppression can 
be extracted directly and possible constraints on QDD effec- 
tiveness can be accurately identified. Each single-axis error is 
obtained from the evolution operator, Uq' Ni \ by projecting 
along the particular axis of interest and performing a partial 
trace over the system. The order of error suppression can then 
be quantified by the scaling of the single-axis error as a func- 
tion of either the total evolution time or the minimum pulse 
interval. We choose the minimum pulse interval since experi- 
mentally this quantity is always lower-bounded, and plays an 
important role in the ultimate performance limits of UDD OTI 
and DD in general l32l . 



The resulting QDD evolution operator, U i 



(JV 2 ,iVi) 

n > 



contains 



all the information regarding decoherence suppression for 
each single-axis error. Construction of the final evolution op- 
erator is accomplished by first considering the inner sequence 
evolution, Z (n ^(t). Let Eq. ^ be partitioned such that 
H = H + + H-, where 



B x 



B,, 



and 



H- = I(g)Bi + tj z (£)B z 



(11) 



(12) 



Clearly, the element of the MOOS comprising Z^ Nl '(r) anti- 
commutes with H + and commutes with #_, i.e., [H±, Z]± = 
0, where the plus and minus sign subscripts signify the anti- 
commutator and commutator, respectively. 

The Z-type UDD sequence Z^ Ni '{t) is effective against 
the anti-commuting Hamiltonian, H + . Z^ Ni \t) is com- 
pletely ineffective against unwanted interactions within H- . 
Any additional errors associated with H- must be addressed 
using another member of the MOOS. The inner sequence evo- 
lution can be expanded in terms of H± ll24ll . 

. . . e -i[-H++H-]s2T e -i[H + +H-]a 1 T ^) 

where the anti-commuting and commuting properties of H±, 
respectively, have been used. Transforming into the interac- 
tion picture with respect to H-, we can write Z^ Ni \t) — 



U 



{t)U { z Ni) (t), such that f/i JVi; (r) = e 



(JVi), 



T and 



f«p(-*y f z (t)H+(t)dt\ . (H) 



The modulation function f z {t) = (— l) 3 1 is defined for 

^ ^ Ei=i s ^ r ' X^=i s i T \ an( l i I s tne time-ordering oper- 
ator. H + in the rotating frame with respect to i?_ takes on the 
form of a power series expansion in t, 



(15) 



k=0 



The power series form of H + (t) is useful (though not essential 
1271 ) for the proof of UDD and therefore the suppression of 
error associated with H+ I13[3TJ|33]|34)- All constants of the 

(k) 

expansion are condensed within H\_ , along with the fc-fold 
commutator 

[ k H_,H+] = [H-,[H-,--.[H-,H+]--]]. (16) 

Using time-dependent perturbation theory, ui Nl \r) is ex- 
panded in the Dyson series 

oo 

(t) = J2Y, H + n) ■ ■ ■ H { + kl) F^ (r), (17) 

n=0 k„ 
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where k n = {Aii, k n } with fej = 0, 1, ... for all i, and all of 
the time-dependence of the expansion has been placed in 

F^(r) = (-i) n / •••/ If <//;/•:/;< 

Jo Jo Jo j =1 

(18) 

The proof of UDD is completed by parametrizing tj as tj — 
rsin 2 (0j/2), Fourier expanding f z (tj), and showing that 

Fz ' k "' > (r) = for all odd values of n when n + 53?= i &j — 
TVi E4l . All even orders of the expansion are proportional to 
unity or a z , and are therefore not associated with H + . The 
expansion ultimately yields 

Z^ Ni \t) = e -^-MS <JVl) r+0((r||//;(r)||) Jv i + 1 ) ; (19) 

where H' + (r) is a generic single-qubit system-bath Hamilto- 
nian and 

H'_ (t) = B'j (t) +<t z ®B' z (t) (20) 

is composed of environment operators dependent on the min- 
imum pulse interval. Note that these environment operators 
are not the same operators defined in Eq. Q, but combina- 
tions of the original Hamiltonian operators resulting from the 
perturbation expansion. 

Left with only dephasing errors, terms proportional to a z , 
the process is continued again by defining H + = a z <£> B' z {t) 
and H- = B'j(t) for the outer X-type UDD sequence of 
Eq. d5Jl. H± is defined in this way such that [H±,X]± = 
0, a similar condition to that required for the inner Z-type 
sequence. The resulting evolution of Eq. <|5j can be summed 
up as 

u (N lt N 2 ) =B ,, ir)+ £ a^B'^r). (21) 

H<£{x,y,z} 

Once again, the time-dependent environment operators, 
B'^{t), differ from the previously defined environment opera- 
tors. Each B'^(t), /i £ {x, y, z}, contains the uncompensated 
decoherence along each of the qubit axes. Bounds on the or- 
der of error suppression are derived analytically in Ref. (35). 



the sum of singular values of A. (The choice of norm is some- 
what arbitrary; we could have used any other unitarily invari- 
ant norm ||36l .) Thus EAt) > 0, and we are interested in 
how the single-axis errors scale as a function of the minimum 
pulse interval. 



C. Distance measure 

The advantage of the single-axis error over other measures 
such as polarization or fidelity is that these typically scale 
with the overall minimum order of decoherence suppression 
and therefore do not provide detailed information about the 
structure of the unitary evolution operator itself. However, 
the single-axis errors of course do not tell the whole story of 
QDD performance. A useful overall fidelity-loss measure is 
the distance [37] 

D(U,G) = ^_L= r>am\\U-G®§\\ F , (25) 
y/d s d B * 

Here ds and ds denote the dimensions of the system and bath 
Hilbert spaces, respectively, U is the actual system-bath uni- 
tary evolution operator, G is the desired system-only unitary 
operator, and $ is a bath operator. Grace et al. l37l give an 
explicit form for this distance measure, so that in our numeri- 
cal simulations we do not need to compute the minimum over 
$. In our case G = I is the desired system unitary operator 
since the goal of DD is to remove the system-environment in- 
teraction while effectively acting trivially on the system. The 
advantage of the distance measure D(U, G) of Eq. ([25} over 
the standard Uhlman fidelity or trace-norm distance [38 1 is 
that it is state-independent and can be correlated directly with 
the results obtained for the single-axis errors. 

We shall show that D(U,I) scales in the same manner as 
min^ E/j,(t). It will also illuminate some interesting features 
of QDD based on the inner sequence order not captured by the 
single-axis errors. 



D. Scaling 



B. Single-axis errors 

One of our goals is to characterize the performance of QDD 
with respect to the remaining system-environment interaction 
operators, B"(r). The error is quantified by what we refer to 
as the single-axis error £„: 



E,(t) = \\B"(t)\\ f , 



(22) 



where 



B' 



^(T)=Tr s (^ % V), (23) 
and where is the Frobenius norm of A, i.e., 

\\A\\ F = TyVa^A, (24) 



We parametrize the strength of the pure environment dy- 
namics and system-environment interactions, respectively, as 



\H F 



J 



SB | 



Jot 



\B n 



(26) 



where ||A|| is the standard sup-operator norm, namely the 
largest singular value (largest eigenvalue of V A^A): 



\\A\\ =sup 



w \/(V#> 



(27) 



The effectiveness of DD tends to be greater in the regime 
where the environment is essentially static and the duration 
of the free evolution is much smaller than the environment 
correlation time: Jr < 1 and /3 <C J. 
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We model the environment as a four-qubit bath with opera- 
tors 



with 



(28) 



characterizing the dynamics of the bath and system-bath in- 
teractions. The operators are composed of one- and 
two-body terms, where i, j index the bath qubits, /i, a, € 
{1,1, y, }, where a 1 = I is the 2x2 identity matrix, and 
c a/3 e 1] are coefficients chosen uniformly at random. 
Constructing £> M in this manner permits general two- and 
three-body interactions between the system and the environ- 
ment, and also facilitates a direct comparison with Ref. 11261 . 
where the same model was used. 

The single-axis error (t) will be dominated by the low- 
est non-vanishing order of r, which we denote by n M . I.e., 
En ~ 0(r n ») provided the first — 1 terms of the power 
series expansion of B'^(t) vanish. More expliclity, we write 



K(r) - E K 



where 



"U) - \' , 



.fx r> r> r> 



(29) 



(30) 



and where dj = {cci, ay} such that ctj € {l,x,y,z} 
and \i = Y[j a j> with the Pauli product rules xy — z, zx — 
y, yz = x. In this manner we have a convenient notation to 
identify the environment operators present in each single-axis 
error. For example, if E z (t) ~ 0(t 2 ) then the possible sum- 

mands comprising J?"^ 2 ' will be proportional to BqB z , B z Bq, 
B x B y , and B y B x . The constituent operators, B a ., are the 
environment operators initially defined in Eq. (j4j). The only 
terms of interest are those proportional to t™'*, since these 
are the dominant terms of B'^(t). Using Eqs. ( |22| i, ( |29"| > and 
( 30 1, and the submultiplicativity property of unitafily invariant 
norms [36], the single-axis error is 

E^r) ~ \\B'^T n »\\ F 

< r n ^\r^ a J\\B ai \\ F ---\\B an J F 

= r n ^<-a n M---J a ^, (3D 



where we have only kept the leading order contribution in 
t. The coupling strength parameters defined in Eq. (26 1 
have been incorporated into the sum, such that J\ = (3 and 
Jx,y.z < J- The parameters r£, a a account for a con- 
version factor between the Frobenius and sup-operator norms. 
Note that the reason we chose to work with the Frobenius 



norm is that the distance measure (25 i is expressed in terms 
of this norm. Factoring J"^ out from the sum, the desired 
functional form of the single-axis error is 



(33) 



and with 



la = Ja/J < 1. 



(34) 



We have written the single-axis error as Eq. p2\ in anticipa- 
tion of our numerical results, where we plot En as a function 
of the dimensionless parameter Jr. In the /3 <C J regime it is 
J which sets the relevant bath timesale and hence we expect 
that Jt < 1 should be a necessary condition for DD to be 
beneficial over uncontrolled free evolution, and we shall see 
that our simulations support this expectation. The quantity \n 
does not depend on r and hence will play the role of a constant 
offset. In the next section we shall unravel the connection be- 
tween the suppression order n„ and the sequence orders N± 
and N2 ■ 



V. NUMERICAL RESULTS 

We now present a numerical analysis of QDD based on the 
single-axis errors and the overall distance measure D(U, G). 
The initial focus is the single-axis error, which we use to quan- 
tify the decoherence suppression as a function of inner and 
outer sequence orders. 



A. Single-axis errors 

Figures [T] and [2] display the single-axis errors as a function 
of Jt for 7V 2 = 3 and N 2 = 4, respectively, with the inner 
sequence order varying from Ni = 1 to Ni = 10. The pure 
bath and system-bath interaction strengths were adjusted such 
that the Hamiltonian is dominated by the interaction, /3 <C J. 
Each data point corresponds to a single cycle of the QDD se- 
quence averaged over 50 random instances of the parameters 



log(^) - n,, log(Jr) + log . 



(32) 



c a p appearing in Eq. (28 1. Since we keep the minimum pulse 
interval fixed, the totaTsequence duration increases with in- 
creasing N\ and The first thing to notice about Figures ([T]) 
and Q is that they match the prediction of Eq. ( 32 1 very well 
in the regime of small Jr. Namely, in all cases we observe a 
constant slope, until Jt ~ 0(1). This is also in agreement 
with the result of Ref. Il26l . 

A summary of the scalings for E^Jt) are given in Table [I] 
for all combinations of Ni,N 2 € {1,2,..., 10}. The val- 
ues of were extracted by performing linear regressions for 
log 10 [Eu(Jt)] between log 10 (Jt) = —9 and the values of 
log 10 (Jr) indicated by the vertical lines in Figs. [T] and [2] and 
rounding to the nearest integer (in all cases the deviation from 
an integer value was at most in the third significant digit). We 
shall return to Table[l]after presenting and discussing the data 
in the figures. 

Let us then consider in detail the effect of varying the in- 
ner and outer sequence orders on the single-axis errors. When 
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Log,„(Jr) Log 1( ,(jT) 

FIG. 1: (color online) Single-axis errors after one cycle of l N 
for outer sequence order iV2 = 3 and inner sequence orders iVi = 
1, 2, . . . , 10 as a function of Jr, averaged over 50 random real- 
izations of the bath operators _B M . Error bars are shown but are 
very small. In all our simulations we set J = 10~ 4 and /3 = 
10~ 6 . Single -axis error values were computed for log 10 (Jr) = 
—9, —8, . . . , 2. Lines are guides to the eye. E x is the solid green 
line, E y is the dotted red line, and E z is the dashed black line. Note 
that log(i? 2 ) is the same in all six plots, with a slope of N2 + 1. The 
slope of \og(E x ), on the other hand, is N\ + 1. The slope of log(E y ) 
is N 1 + 2. Vertical lines denote the largest value of Jr utilized in the 
linear regression used to extract the slope ro M . 



N 2 > Ni, higher order suppression is expected for the er- 
rors that correspond to the system basis operators which anti- 
commute with the member of the MOOS comprising the outer 
X-type sequence. Thus the single-axis errors E v (t) and 
E z (t) are most heavily suppressed. 

Since only the outer sequence can suppress z-axis, or Z- 
type errors [recall Eq. (|5j], E z (Jt) only gains additional er- 
ror suppression if the outer sequence order is increased. In 
Fig. Q, N 2 = 3 and E z (Jt) ~ <D[( Jr) 4 ] for all Ni, exhibit- 
ing error suppression of the first N 2 terms of the interactions 
proportional to a z . Thus QDD operates with UDD efficiency 




LogioGM Log 10 (Jr) 

FIG. 2: (color online) Single-axis errors after one cycle of U^ 1 '^ 
for outer sequence order N2 = 4: and inner sequence orders Ni = 
1, 2, . . . , 10 (left to right, top to bottom) as a function of Jr, aver- 
aged over 50 random realizations of the bath operators B M . Other 
details as in Fig. [TJ except that the single-axis error suppressed by 
both the inner and outer sequence, E y (r), exhibits a strong depen- 
dence on the parity of the inner sequence. 



for error suppression by the outer nested sequence alone. 

In a similar manner to E z (Jt), the behavior of E x (Jt) 
can be attributed to one of the two nested sequences. Namely, 
the error measured by E x is associated with a x , which only 
anti-commutes with the MOOS operator present in the inner 
Z-type sequence. Determined solely by the inner sequence or- 
der, E x (Jt) ~ 0[(Jt) 1+ ]. Essentially, the outer sequence 
has no effect on the order of error suppression for E x (Jt), as 
can be seen from Fig.[T]for iV 2 = 3. 

The interpretation for E y (Jr) is not as simple, since this 
single-axis error is compensated by both the inner and outer 
sequences. One might expect both the inner and outer se- 
quence to contribute to F-type error suppression, i.e., E y (Jr) 
to scale with ( j T ) m ax(Ar 1 ,Af 2 )+i^ fi owever; if t hi s W ere the 
case then, e.g., the case N\ = 1 would display an equal or- 
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TABLE I: Summary of the scaling for all single-axis errors. Values of were extracted by performing a linear regression, rounded to 
the nearest integer, fitting the slopes of the straight line portions of the curves displayed in Figs. [T] and [2] and the additional Figs. |7|10] in 
Appendix|X] between log 10 (Jr) = —9 and the values of log 10 ( Jr) indicated by the vertical lines in these figures, (a) E x , (b) E y , and (c) 
E z for Ni,N 2 £ {1, 2, . . . , 10}. For n y and n z the outer sequence order N 2 is displayed in the top row and the inner sequence order Ni in 
the first column. Each of the single-axis errors is dominated by the lowest order of Jr, denoted n M , therefore E^ ~ 0[(Jr)" f ']. Additional 
simulations (not shown) fully continue the trends seen in this table and summarized in Eqs. |35}-l[37} all the way up to Ni,N 2 < 24. 



der of error suppression for both E y (r) and E z (t). Instead 
we find that E y (Jr) ~ 0[( Jt) Ni+2 } for N 2 " = 3. Thus 
the suppression of E y (Jr) is constrained by Ni, even when 
N2 > N±, though it is larger by one order of magnitude than 
UDD error suppression efficiency for the inner sequence. 

Similar observations apply for all odd-order outer se- 
quences we have analyzed (see Appendix [A] Figures [7] [9 



Comparing the case of the anti-symmetric outer sequence 
of Fig. [T] to that of the symmetric sequence of A^ — A in 
Fig. [2] one notices immediately that there is a qualitative dif- 
ference. The single-axis error E y (Jr), the component anti- 
commuting with both the inner and outer sequences, a v , fluc- 
tuates strongly as a function of Ni. A similar effect is ob- 
served for other even values of (see Appendix [A] Fig- 



TT] and \m. Odd-order sequences are anti-symmetric with ures [TO] [12| and [R). Only the outer sequence order has 



respect to time-reversal, and the conclusions concerning the 
case N2 — 3 can be generalized as follows: when the outer 
sequence is anti-symmetric, terms in the QDD evolution oper- 
ator which anti -commute with only one element of the MOOS 
are suppressed with UDD efficiency, determined by the order 
of the nested sequence composed of the corresponding anti- 
commuting MOOS operator (this applies to a x and er 2 ). In 
contrast, terma that anti-commutes with both elements of the 
MOOS are suppressed to one order beyond UDD efficiency, 
dictated exclusively by the inner sequence order (this applies 
to a y ). Below we will see how this observation is modified 
when we consider larger values of A^. 



been changed, therefore this characteristic is entirely depen- 
dent on the fact that the outer sequence is now symmetric. 

Analogous to the anti-symmetric outer sequence, Fig. [2] 
shows that the single-axis error E z (Jt) is independent of 
the inner sequence order. The scaling E z (Jt) ~ C[(Jr) 5 ] 
holds for all Thus the single-axis error E z (Jt) again ex- 
hibits UDD efficiency, independent of the parity of the outer 
sequence. Similarly, again E x (Jt) ~ 0[{ Jt) Ni+1 ] in Fig. [2] 
However, when we consider the n z results for all values of 
Ni and N2 we find that there are exceptions to this simple 
behavior. As can be seen from Table |IJ when = 1 and 
N2 > 4, the value of n z is fixed at 4. The same phenomenon 
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is observed for Ni = 3 and N 2 > 8. 

On the basis of our numerical data we can summarize the 
scaling of the X and Z-type single-axis errors as follows: 



and 



iV 2 + l 
N 2 + 1 
27Vi +2 



Nx even 

iVi odd, N 2 < 2N X 
Ni odd, N 2 > 2N X 



(35) 



(36) 



Qualitatively, we expect that when the inner sequence 
works imperfectly, as is the case for N\ odd, the lowest order 
sequence will determine the scaling of the single axis error, 
and this is what is stated in Eq. d36j. 

As is clear from Fig. (pjl, the scaling of E y (Jr) is de- 
pendent on the parity of Ni. If the inner sequence is of 
odd parity then E y (Jr) ~ 0[{ Jt) Ni+2 ] when N 2 is odd as 
well, or Ey(Jr) ~ 0[{ Jt) Ni+1 ] when N 2 is even. Thus 
the scaling of E y is dominated by the inner sequence or- 
der when Ni is odd. The situation changes when Ni is 
even. Now, if N 2 is odd the sequence is still anti-symmetric, 
however there is an immediate improvement in error sup- 
pression, E v {Jt) ~ 0[(j T )max(JV 1+ i,A' 2 )+i^ if t h e com- 
plete sequence is fully symmetric (both N\ and N 2 even) we 
also find a scaling dependent on both the inner and outer se- 
quence orders, E y (Jr) ~ 0[( j T )max(JV 1 ,JV 2 )+ij_ We thus 

see that the suppression of interactions which anti-commute 
with both elements of the MOOS depends sensitively on the 
parity of the inner sequence order, and is summarized for 

Ey(jT) ~ 0[{Jr) n y] as 



max(7V 1 ,7V 2 ) + 1 
max(iVi + 1,N 2 ) + 1 
N t + 1 
Ni +2 



Ni even, N 2 even 
N\ even, N 2 odd 
Ni odd, N 2 even 
TVi odd, N 2 odd 



The dependence upon the symmetry of the inner sequence, 
the parity of Ni, was first noted by Wang & Liu in the context 
of overall QDD performance [27). The dependence on the 
outer sequence symmetry, however, was not noted previously. 
Our results show that the symmetry of the outer sequence im- 
pacts the efficiency of a y error suppression as well. 

The efficiency of QDD error suppression is directly related 
to the efficiency of UDD. Each interaction that anti-commutes 
with at least one member of the MOOS is expected to achieve 
UDD efficiency. Fully symmetric QDD, i.e., even order N% 
and N 2 , recovers the efficiency of UDD for all single-axis er- 
rors. Consequently, QDD performs with optimal efficiency 
when it is fully symmetric. 

The interactions addressed only by the inner or outer se- 
quence are separately suppressed with UDD efficiency corre- 
sponding to the order of the corresponding sequence perform- 
ing the decoherence suppression. Equivalently, interactions 
which anti-commute with only one member of the MOOS 
are suppressed with UDD efficiency in the QDD scheme. 
On the other hand, error suppression of the interactions anti- 
commuting with both elements of the MOOS is dependent 
upon the parity of both the inner and outer sequence orders. 



B. Intermediate single-axis errors 

Rather than consider the single-axis errors just at the end 
of the QDD sequence, here we consider the single-axis er- 
rors prior to the application of each X-type outer sequence 
pulse. We will refer to these as "intermediate single-axis er- 
rors" since they are extracted during the QDD evolution, un- 
like those presented in Figures [T] [2] and [7p4] which are ex- 
tracted at the end of the complete evolution, y studying this 
intermediate time-dependence of the errors we shall gain an- 
other interesting perspective on the manner in which the QDD 
sequence suppresses decoherence. 

Let us define a set of "intermediate QDD" sequences as 

jjgfij) _ Z ^\a j T)XZ^(a S - 1 T) ■ ■■XZ^\s 1 r), 

(38) 

where j € {1, . . . , N 2 + 1}. Thus, except for j = 1, U^ 1 ^ 
contains j — 1 X-type pulses sandwiched between j Z-type 
UDD sequences. When j = 1 



is just the UDD sequence. We also separately define 



U, 



(N u N 2 +2) ^ yiNuNz) _ y N2 f T (N u N 2 +l) 



= x ,y - 2 u' 



(39) 



(40) 



i.e., the complete QDD sequence, Eq. ([5j. Note that 
contains a final X pulse if N 2 is odd, but not if 



(JVi,JV 2 +2) 



N 2 is even. Similarly to the error expansion pi) , we have the 
intermediate error expansion 



B<P(t)+ J2 

p,£{x,y,z} 



^ <8> B®(t). (41) 



In analogy to Eq. ( |22] > we can now define the intermediate 
single-axis errors as 



E<g\r) = \\B^{t)\\ 



(42) 



where j e {l,N 2 + 2}. Note that for odd N 2 the er- 

differ by a single instantaneous 



rors E, 



(JV2+1) 



and E, 



(N 2 +2) 



X pulse (which is significant, as our simulations results will 



(JV 2 + 1) _ £,(IV 2 +2) 



so that 



demonstrate), while for N 2 even E t 
below we do not plot El N2+2 ^ in the even case. 

Figures [3] and |4] display the intermediate single-axis errors 
for N 2 = 3 and N 2 — 4, respectively, with Ni = 1, 2, . . . , 6. 
The coupling parameters are fixed at J = 10~ 4 and /? = 10~ 6 
as in the previous figures. Additional results for odd N 2 are 
given in Appendix[A]in Figuresjl9]and 21 and for even N 2 in 
Appendix [A| in Figures [20] and [22] 

Several features are noteworthy in these figures. 



(i) Ex and E y y ' are equal and substantially smaller than 
E^\ and the difference grows as N\ is increased. This is 
because the inner Z-type sequence only suppresses the X and 
F-type errors, and the point j ' = 1 does not include the first 
X-type outer sequence pulse. Formally, this is expressed by 

E$Ht) = \\BW\\ F ~( ai T) N * +1 , ve{x,y} 
E^(t) = ||B,W|| F ~1, ve{I,z}. (43) 
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FIG. 3: Intermediate single-axis errors for N2 = 3, as defined in 
Eq. l|42}. For given j, the intermediate single-axis error is computed 
after j innner Z-type UDD sequences separated by j — 1 pulses. 
There are A^2 + 1 inner UDD sequences. The point j = 5 is the 
last X pulse at the end of last inner UDD sequence, as required for 
odd N2 • Note that because the data points labeled j = 4 and j = 5 
are separated by a single X pulse, and our pulses are instantaneous, 
these points have no actual time delay between them. 



(ii) The intermediate single-axis errors all fluctuate 
throughout the QDD evolution. This is due to a reshuffling of 
the errors after each outer sequence X-type pulse is applied, 
a simple consequence of the rules of Pauli matrix multiplica- 
tion. To see why in some detail, consider the effect of the first 
X = a x ® / pulse: 



XU. 



n 



B 



(i) 



E 

i>e{x,y,z} 
(1) 



B (}) 



(44) 



B 



a" 



B 



(i) 



B 



(i) 



where we dropped factors of i. The reshuffling effect is clear: 
for example, the error single-axis z-type error now comes 
from By 1 ^. To explain the j = 2 behavior we should consider 
the effect of multiplying XUq' 1 ' by the next inner UDD 
sequence Z^ Nl ^> (s 2 t). The ith inner UDD sequence has the 
expansion 

Z( Nl \ Si T) = B Iti +a x ®B Xti +a y ®By t i+a z ®B Zti , (45) 
where similarly to Eq. ( |43| we have 
115,, 



■>H,i\\F 
\BvittF 



1, ue{I,z}. (46) 



Using this to carry out the multiplication to the next order we 
have 



^,2) = z(Nl) 



= [Bl,2 



■f CT X <E) B x>2 + ( 

+ g x ® M 1} + < 



(47) 

B z ,2] 



Consequently 



Ei 2 \r) = \\B l!2 B^+B x , 2 BP+B y , 2 B^+B z , 2 BP\\ F 



E y 2 \r) = \\B Ia B^+B x ^+B ya B^+B za B { p\\ F 
- 1 

Ei 2 \r) = \\B I>2 BW + B X>2 B^ + By . 2 B ( p + B Z . 2 B^\\ F 
-2( Sl r) Ari+1 +2( S2 T) JVl + 1 , (48) 



where (r) is dominated by Bj^Bj 1 ' and E y z> (r) is dom- 
inated by Bj t2 Bz, neither of which is suppressed, whence 
the ~ 1 result. On the other hand every one of the terms in 
E^ (r) is suppressed. Hence, as can be seen in Figuresland 
|4] (and their companions, Figures 19p2 in the Appendix), at 
j = 2 both the X and y-type errors have increased relative to 
j = 1, while the Z-type error has decreased. One can simi- 
larly understand the remaining oscillations of the intermediate 
single-axis errors in terms of this reshuffling of error types. 

(iii) E x and E z oscillate out of phase, while E y oscillates 
in phase with E x for even N 2 , but not necessarily for odd N\. 
This is again a consequence of error reshuffling. The F-type 
error behaves differently from the other two since it experi- 
ences suppression from both the inner and outer sequences. 
For the same reason we always find E y < E x . 

(iv) E x attains its minimum for j = 1 and then slowly in- 
creases, though while maintaining its suppression order. This 
is because the X-type error is suppressed only by the inner 
sequences, and these are simply applied to it with fixed order 
(A?i), a total of N 2 or N 2 + 1 times. Repeated application of 
the inner UDD sequence is similar to the periodic DD (PDD) 
protocol, whose performance is well known to deteriorate as 
time grows lfl2l[39l . The reason is that the error accumulates 
over time, without a mechanism for reducing it. 

(v) There does not appear to be much of a difference be- 
tween even and odd values of N 2 in terms of the intermediate 
single-axis errors. One difference is that E y ' tends to be more 
erratic for odd N 2 at high j values. We do not have a simple 
explanation for this behavior. Another difference is that for 
even N 2 all single-axis errors have the same final value when 
Ni = N 2 , but for odd N 2 the X-type error is always slightly 
worse at the end of the sequence, thus setting the bottleneck. 
Perhaps additional pulse interval optimization can remove this 
asymmetry. 

(vi) Only at the very end are all three single-axis errors si- 
multaneously small. Thus, while suppression of one error type 
can be achieved in the middle of the QDD sequence, one must 
wait until its completion to suppress all errors. 
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TABLE II: Summary of the scaling of the overall distance measure D with respect to inner and outer sequence orders, Ni and A 2 , respectively. 
Values of nn were extracted by performing a linear regression, rounded to the nearest integer, fitting the slopes of the straight line portions 
of the curves displayed in Figs. [T] and [2] and the additional Figs. 15|18 in Appendix |a| between log(JY) = —9 and the values of log(Jr) 
indicated by the vertical lines in these figures. The outer sequence order N 2 is displayed in the top row and the inner sequence order TVi in 
the first column. We find that, as expected, no = min(na;, n y , n z ). Additional simulations (not shown) fully continue the trends seen in this 
table and summarized in Eq. \50\ all the way up to Ni,N 2 < 24. 
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FIG. 4: (color online) Intermediate single-axis errors for = 4. 
As in Fig.[3]except that there is no final X pulse for N 2 even, i.e., 
there are N 2 + 1 inner UDD sequences separated by X pulses. 



C. Overall performance 



While the single-axis error analysis presented in the previ- 
uous two subsections helps in unravelling the mechanism of 
QDD performance, it does of course not tell the whole story. 
We now present our results for the distance measure D(U, I) 
[Eq. (|25j>], which provides a complete quantitative description 
of QDD performance. We expect this overall performance of 
QDD to be dictated by the lowest order of r present in the 
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FIG. 5: (color online) Overall QDD distance measure after one cycle 
of Ufo 1 ' for outer sequence order A2 = 3 and inner sequence or- 
ders Ni = 1, 2, . . . , 6, as a function of Jr, averaged over 50 random 
realizations of the bath operators B^. The performance of QDD pro- 
gressively improves with increasing Ni up to N 2 — Ni, indicating 
that min{Ai, N 2 } dominates QDD performance. 



final evolution operator, Eq. pT[ ), i.e., 

D - C[(Jr)" D ] 

where 

n D = mm(n x ,n y ,n z 



(49) 



(50) 



Overall QDD performance for = 3 and N2 — 4 is 
shown in Figures [5] and |6] respectively. The outer sequence 
order N2 is fixed and the inner sequence order Ni is varied 
from 1 to 6. These results are for the same model considered 
in the previous subsection. Additional results are given in Ap- 
pendix|A]for N 2 — 1, 2, 5, 6 (see Figures T3JT8 1. A summary 
of the distance scaling results is presented in Table [II] 
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FIG. 6: (color online) Overall QDD distance measure after one cycle 
of Uq 1,N for outer sequence order N2 = 4 and inner sequence 
orders Ni = 1, 2, . . . , 6, as a function of Jt, averaged over 50 ran- 
dom realizations of the bath operators B^. The dependence of the 
order of error suppression on xam{N\, N2} is again observed. 



Considering iV 2 = 3 first (Fig. E), when A^ 2 > Ni the 
overall order of error suppression is hindered by the inner 
sequence order. This is evident by the increasing order of 
error suppression as N\ increases. In this regime the lower 
sequence order is determined by the inner sequence, there- 
fore the scaling of D is equivalent to that of E x , i.e., D ~ 
0[( Jt) Ni+1 ]. As iVi passes N2, there is a saturation of error 
suppression corresponding to a performance bounded by the 
lower outer sequence order. The amplitude of performance 
increases slightly beyond jVjj = Ni, however begins to de- 
crease when Ni > N2, as evidenced not by the slope but by 
the offset of the distance curves. Namely, the ordering, from 
worst to best, is Ni — 6, 5, 4. The latter is an interesting fea- 
ture not easily deduced from the single-axis errors. Increasing 
the inner sequence order results in an accumulation of error 
for the single-axis error dominating the performance; when 
Ni > N2 + 1 this corresponds to E z . 

The results are similar for N2 — 4, as shown by Fig. [6] 
The order of error suppression, given by the slope, increases 
until N\ = N2 in correspondence with an overall performance 
dominated by the lowest order of r present in E^. In addition 
to the saturation of the order of error suppression, we again 
observe an offset-related deterioration. Namely, Ni = 6 is 
slightly worse than Ni = 5. 



VI. CONCLUSIONS 



axis errors in the regime where DD is most beneficial. We 
constructed our QDD sequences with N\ Z-type pulses com- 
prising the inner sequence, and N 2 X -type pulses comprising 
the outer sequence. We found that the system-bath interaction 
term proportional to a x is suppressed with UDD efficiency for 
all values of N\ and 7V 2 [Eq. ( [35] )] . The interactions propor- 
tional to a z and <j v both exhibit parity effects [Eqs. ( 36 1, ( 37 1] 
whose origins are the symmetry or anti-symmetry of the inner 
and outer UDD sequences. Of course, permuting the pulse 
types of the inner and outer sequences will correspondingly 
modify these conclusions. 

We also performed an analysis of the intermediate time- 
dependent performance of QDD. We found that the single- 
axis errors are strongly time-dependent, oscillating between 
outer-sequence pulses, until they all converge to nearly the 
same value after the final outer-sequence pulse. The closest 
convergence occurs for QDD sequences with equal inner and 
outer orders. 

Finally, we computed the overall performance of QDD us- 
ing an appropriate distance measure, and reconciled its scal- 
ing with that of the single-axis errors. We showed that overall 
QDD error suppression scales with the lowest order of single- 
axis error suppression, i.e., the first non-vanishing contribu- 
tion appears at order min(iVi,iV 2 ) + 1. QDD accomplishes 
this by applying (Ni + l)(A r 2 + 1) pulses. We conjecture that 
similarly, for NUDD with K nested UDD sequences using 
(Ni + l)(N 2 + 1) • • • (N K + 1) pulses, the first non-vanishing 
contribution will appear at order mm.j(Nj) + 1. 

In this work we treated the pulses as ideal, instantaneous 
operations. However, this is of course an idealization. An 
important topic for future study is robustness with respect to 
pulse errors, whether random or systematic. This topic has 
been addressed for UDD both theoretically B0ll4T1 and exper- 
imentally l42l . and the overall conclusion is that pulse errors 
can have a dramatic negative impact unless they are compen- 
sated for. Some combination of pulse shaping and optimiza- 
tion will surely be required to overcome this problem in the 
context of QDD as well. 



Acknowledgments 

We are grateful to Gonzalo Alvarez, Wan-Jung Kuo, Ste- 
fano Pasini, Dieter Suter, and Gotz Uhrig for very helpful 
discussions. DAL acknowledges support from the U.S. De- 
partment of Defense and the NSF under Grants No. CHM- 
1037992 and CHM-924318. 



This work presents a comprehensive numerical analysis 
of the error suppression characteristics of QDD. This was 
achieved by isolating the single-axis errors associated with 
each system basis operator in the system-bath interaction. The 
order of error suppression was determined by computing the 
single-axis error as a function of the minimum pulse inter- 
val. We performed our analysis for a model in which the 
system-environment interaction dominated the internal bath 
dynamics, so that we could study the properties of the single- 
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Appendix A: Additional Numerical Results 



In this appendix we present additional figures in support of 
the numerical results presented in the body of the paper. 
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FIG. 7: (color online) Single-axis errors after one cycle for iV2 = 1 
and A^i = 1, 2, . . . , 10 (left to right, top to bottom). See Fig.[T]for 
additional details. 




FIG. 8: (color online) Single-axis errors after one cycle for N2 — 2 FIG. 9: (color online) Single-axis errors after one cycle for N2 — 5 
and Ni = 1, 2, . . . , 10 (left to right, top to bottom). See Fig.|2]for and Ni = 1, 2, . . . , 10 (left to right, top to bottom). See Fig.[T]for 
additional details. additional details. 
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FIG. 10: (color online) Single-axis errors after one cycle for N2 = 6 FIG. 1 1 : (color online) Single-axis errors after one cycle for N2 = 7 
and Ni = 1, 2, . . . , 10 (left to right, top to bottom). See Fig.|2]for and Ni = 1, 2, . . . , 10 (left to right, top to bottom). See Fig.[2]for 
additional details. additional details. 
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FIG. 12: (color online) Single-axis errors after one cycle for N2 = 8 FIG. 13: (color online) Single-axis errors after one cycle for N2 = 9 
and Ni = 1, 2, . . . , 10 (left to right, top to bottom). See Fig.|2]for and Ni = 1, 2, . . . , 10 (left to right, top to bottom). See Fig.[2]for 
additional details. additional details. 



16 




-8 -6 -4 -2 2 -8 -6 -4 -2 2 

Log,„(Jr) Log 1( ,(jT) 



FIG. 14: (color online) Single-axis errors after one cycle for N2 = 
10 and Ni = 1, 2, . . . , 10 (left to right, top to bottom). See Fig. [2] 
for additional details. 
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FIG. 16: (color online) Overall QDD distance measure after one cy- 
cle for N 2 = 2 and Ni = 1, 2, . . . , 6 (left to right, top to bottom). 
See Fig.|6]for additional details. 
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FIG. 17: (color online) Overall QDD distance measure after one cy- 
cle for N2 = 5 and Ni = 1, 2, . . . , 6 (left to right, top to bottom). 
See Fig.|5]for additional details. 
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FIG. 15: (color online) Overall QDD distance measure after one cy- 
cle for N 2 = 1 and Ni = 1, 2, . . . , 6 (left to right, top to bottom). 
See Fig.[5]for additional details. 



FIG. 18: (color online) Overall QDD distance measure after one cy- 
cle for N 2 = 6 and Ni = 1, 2, . . . , 6 (left to right, top to bottom). 
See Fig.|6]for additional details. 
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FIG. 19: (color online) Intermediate single-axis errors for N2 = 1. 
See Fig.[3]for additional details 
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FIG. 21: (color online) Intermediate single-axis errors for N2 = 5. 
See Fig.[3]for additional details 
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FIG. 20: (color online) Intermediate single-axis errors for N2 = 2. 
See Fig.[4]for additional details 



FIG. 22: (color online) Intermediate single-axis errors for N2 = 6. 
See Fig.|4]for additional details 



